Nonlocal corrections to the Boltzmann equation for dense Fermi systems 
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A kinetic equation which combines the quasiparticle drift of Landau's equation with a dissipation 
governed by a nonlocal and noninstant scattering integral in the spirit of Snider's equation for 
gases is derived. Consequent balance equations for the density, momentum and energy include 
quasiparticle contributions and the second order quantum virial corrections. The medium effects on 
binary collisions are shown to mediate the latent heat, i.e., an energy conversion between correlation 
and thermal energy. An implementation to heavy ion collisions is discussed. 



I. INTRODUCTION 

One of the goals of the recent nuclear physics is to 
find the equation of state of nuclear matter. Indeed, the 
dependence of the pressure on the density of nucleons is 
a crucial input for a hydrodynamical modeling of heavy 
ion collisions or of astrophysical events like the big bang, 
supernova explosions and neutron stars . 

In the absence of any direct measurement, it is hoped 
that the equation of state can be deduced from heavy ion 
collisions via the following scheme. Heavy ion collision 
data are fitted with the Boltzmann equation (BE) 
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Arguments of distributions / and energies e are shortened 
as fx = f a (k,r,t), f 2 = fb{p,r,t), f 3 = f a {k-q,r,t), and 
fx = fb(p + Q, r ,t), with momenta k,p,q, coordinate r, 
time t, and spin and isospin a, b. Once the differential 
cross sections \T\ 2 and the functional dependence of en- 
ergy £ on the distribution / are fitted, the equation of 
state is evaluated from the kinetic equation. 

This scheme has two drawbacks. First, accessible 
fits of the quasiparticle energy e are not sufficiently 
reliable since two possible fits, momentum-dependent 
and momentum-independent, result in very contradictory 
predictions giving hard and soft equations of state, re- 
spectively j| . When this more or less technical problem 
is resolved in future, one has to face the second draw- 
back: the BE is not thermodynamically consistent with 
virial corrections to the equation of state. This problem 
is principal for "how can one infer the equation of state 
from the BE if the two equations are not consistent?" . A 
consistency between the kinetic and the thermodynamic 
theories is a general question for the quantum statistics 



exceeding the merits of the nuclear matter. Here we ap- 
proach this question from nonequilibrium Green's func- 
tions. It is shown that the consistency is achieved by a 
consistent treatment of the quasiclassical limit which re- 
sults in nonlocal and noninstant corrections to the scat- 
tering integral of the BE. 

The need of nonlocal corrections can be seen on the 
classical gas of hard spheres. In the scattering integral of 
(|l|) , all space arguments of the distributions are identical, 
i.e., colliding particles a and b are at the same space 
point r. In reality, these particles are displaced by the 
sum of their radii. This inconsistency has been noticed 
by Enskog [|| and cured by nonlocal corrections to the 
scattering integral. The equation of state evaluated from 
the kinetic equation with the nonlocal scattering integral 
is of the van der Waals type covering the excluded volume 
HQ. For nuclear matter, Enskog's corrections has been 
first discussed by Malfliet || and recently implemented 
by Kortemayer, Daffm and Bauer ||. 

The noninstant corrections are closer to the chemical 
picture of reacting gases. In the scattering integral of 
(|l|), all time arguments of the distributions are identi- 
cal what implies that the collision is instant. In real- 
ity, the collision has a finite duration which might be 
quite long when two particles form a resonant state. The 
resonant two-particle state behaves as an effective short- 
living molecule. Like in reacting gases g, the presence of 
these molecules reduces the pressure since it reduces the 
number of freely flying particles. The finite duration of 
micleon-nucleon collisions and its thermodynamic conse- 
quences has been for the first time discussed only recently 
by Danielewicz and Pratt jl]]. The noninstant scattering 
integral and its consequencies for the linear response has 
been also discussed for electrons in semiconductors scat- 
tered by resonant levels || . 

Except for dense Fermi systems, the above intuitively 
formulated nonlocal and noninstant corrections has been 
confirmed by systematic approaches. For classical gases, 
this theory was developed already by Bogoliubov and 
Green Obtained gradient contributions to the 

scattering integral are the lowest order terms of the virial 
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expansion in the kinetic equation |Ll[ | . The first quantum 
kinetic equation with nonlocal corrections has been de- 
rived by Snider p^ j. Recently, it has been recognized 
that Snider's equation is not consistent with the second 
order virial corrections to equations of state. A consis- 
tent quantum mechanical theory of the virial corrections 
to the BE has been developed from the multiple scat- 
tering expansion in terms of M0ller operators |l3| and 
confirmed by Balescu's formalism ]l4j . 

Presented treatment extends the nonlocal and non- 
instant corrections to dense Fermi systems. We follow 
Baerwinkel Jl5| ] in starting from nonequilibrium Green's 
functions and keeping all gradient contributions to the 
scattering integral. Baerwinkel's results are limited to 
low densities (to avoid medium effects on binary colli- 
sions) and not consistent (since he uses the quasiparti- 
cle approximation). Here we describe the binary colli- 
sions by the Bcthe-Goldstone T-matrix which includes 
the medium effects. Instead of the quasiparticle approxi- 
mation, the extended quasiparticle approximation is used. 
This extension is sufficient to gain consistency of the ki- 
netic theory with the virial corrections to thermodynamic 
quantities. 



II. EXTENDED QUASIPARTICLE PICTURE 

We start our derivation of the kinetic equation from 
the quasiparticle transport equation first obtained by 
Kadanoff and Baym @0 
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Like in (|l|), quasiparticle distribution /, quasiparticle en- 
ergy e and wave-function renormalization z are functions 
of time t, coordinate r, momentum k and spin and isospin 
a. Self-energy E, taken from nonequilibrium Green's 
function in the notation of Kadanoff and Baym Jt6| , is 
moreover a function of energy u>, however, it enters the 
transport equation only by its value at pole to = e±. 

Particular forms of the quasiparticle energy and the 
scattering integral we derive for a model and an approxi- 
mation used in nuclear matter for heavy ion collisions in 
the non-relativistic energy domain. The system is com- 
posed of protons and neutrons of equal mass m. They 
interact via an instant potential V. We assume no spin- 
flipping mechanism. As common, the self-energy is con- 
structed from the two-particle T-matrix T R in the Bcthc- 
Goldstone approximation as [T^l, 2, 3, 4) = (1- 

S aia2 )T*(l, 2, 3, 4) + -^<5 aiQ2 (T^l, 2, 3, 4)-T fl (l, 2, 4, 3))] 



E<(1,2)=T«(1,3;5,6)T^(7,8;2,4) 
x G > (4,3)G < (5,7)G < (6,8), 



and E> is obtained from (0) by an interchange >^<. 
Here, G's are single-particle Green's functions, numbers 
are cumulative variables, 1 = (ii,7*i,Oi), and bars denote 
internal variables that are integrated over. Before (^) is 
plugged in (Q), it has to be transformed into the mixed 
representation, [off-diagonal elements in spin and isospin 
are excluded, a\ = a-i = a] 
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and all Green's functions in (||), too. 

The self-energy E is a functional of Green's func- 
tions G. This functional E[G] is converted to the func- 
tional of the quasiparticle distribution E £ [/] via the ex- 
tended quasiparticle approximation |i~^ , |l9[ ] [zi = 1 + 
£jReEi„| ] 



^ > 
G^ = ( y ) 2*zi6(w - ei) + Re ^ , 
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where Gi )W = G a (u),k,r,t) and similarly E. Unlike the 
plain quasiparticle approximation (without the second 
term) used by Baerwinkel [^5[ , approximation ([&]) leads 
to the consistent theory. The first term brings the on- 
shell quasiparticle part, the second term is the off-shell 
contribution. 

The off-shell part plays four-fold role. First, it justifies 
the kinetic equation (^]). Equation (j^) has been origi- 
nally derived from the plain quasiparticle approximation 
neglecting the off-shell drift ^ The off-shell part of G < 
in (||) compensates the off-shell drift so that (|f) is re- 
covered without uncontrollable neglects JIt]]. Second, in 
the quasiparticle energy Si = -M— + ReE^ E1 , the off-shell 
part brings contributions that are essential for the cor- 
rect binding energy |||. Third, (|) provides Wi gner s 
distribution p = J^G < asa functional of the quasipar- 
ticle distribution / | ^7| , p0[ . Fourth, in the scattering in- 
tegral of (||) , the off-shell part results in sequential three- 
particle processes with the off-shell propagation between 
the two composing binary processes. Since the three- 
particle processes are beyond the scope of the present 
paper, they are excluded from scattering integral. 



III. NON-LOCAL SCATTERING INTEGRAL 

Now the approximation is specified and we can start 
to simplify the scattering integral. In contrast to previ- 
ous treatments of degenerated systems, we keep all terms 
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linear in gradients. The gradient expansion of the self- 
energy (H|) is a straightforward but tedious task. It re- 
sults in a one nongradient and nineteen gradient terms 
that are analogous to those found within the chemical 
physics p^LO ). All these terms can be recollected into 
a nonlocal and noninstant scattering integral that has 
an intuitively appealing structure of the scattering inte- 
gral in the BE ([!]) with Enskog-type shifts of arguments. fj] 
In agreement with (DMijl, all gradient corrections result 
proportional to derivatives of the scattering phase shift 
<j> = Im hiT^(Sl,k,p,q,t,r), 
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After derivatives, A's are evaluated at the energy shell 
il — > £i + e 2 . The corrected BE with the collected gradi- 
ent terms then reads [A r — j(A 2 + A3 + A4)] 

dfi deidfi _ dei dfi 

dt dk dr dr dk 
dpdq 
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Unlike in (|l]), the subscripts denote shifted arguments: 
/1 = f a (k,r,t), h = h(p,r- A 2 ,t), f 3 = f a (k-q- 
A K ,r-A 3 ,t-A t ), and f 4 = f b (p+q-A K ,r-A 4 ,t-A t ). 

The A's are effective shifts and they represent mean 
values of various nonlocalities of the scattering integral. 
These shifts enter the scattering integral in the form 



^The basic idea of the recollection can be demonstrated on 
the following rearrangement of the gradient approximation of 
a matrix product C(l,2) = A{\, 3)B(3, 2). In the mixed rep- 
resentation C = AB + f - - + ^ffi), 

see |l(j,|l7| . We denote if — | In A and rearrange the 
product as C = A (B + g« - f » - |f M + % bb ). 
The gradient term in brackets can be viewed as a 
linear expansion of B with all arguments shifted as 
C = A{w,k,r,t)B{ u -%k+%,r-%,t+%;). 



known from the theory of gases []3], Ojl4[ ), however, the 
set of shifts is larger due to the medium effects on the bi- 
nary collision that are dominated by the Pauli blocking 
of the internal states of the collision. 

The physical meaning of the A's is best seen on grad- 
ually more complex limiting cases: 

(0) Sending all A's to zero, (0) reduces to the BE ([j]). 

(1) In the classical limit for hard spheres of the diam- 
eter d, the scattering phase shift cj> — * n — \q\d gives 
A4 = A2 = A-d and all other A's are zero. The En- 
skog's nonlocal corrections are thus recovered. 

(ii) For a collision of two isolated particles interacting 
via the single-channel separable potential, the scattering 
phase shift does not depend on r, t, q and k — p, while 
it depends on k + p exclusively via the energy depen- 
dency, ^ c/>(n- ^(k+p) 2 ). Then A £ii f, 2 = and 
A 3i 4 = A t . Since is the center-of-mass veloc- 
ity, the displacements A4 = A3 represent a distance over 
which particles fly together as a molecule. 

(iii) For two isolated particles interacting via a general 
spherical potential, the scattering phase shift reflects the 
translational and the spherical symmetries. From trans- 
lation of the center of mass during A t follows the relation 
for the molecular flight i(A 4 + A 3 - A 2 ) = ^A t . From 
the spherical symmetry follows that the sum of relative 
coordinates of the particles a and b at the end and at the 
beginning of collision has the direction of the transferred 
momentum ^(A 4 — A 3 + A 2 ) = A-d (Enskog-type shift), 
and that the difference has the perpendicular in-plane 
direction ^(A 4 — A3 — A2) = r§E^E|] a (rotation of the 
molecule). The nonlocality of the collision is thus given 
by three scalars At, d and a. 

(iv) For effectively isolated two-particle collisions (no 
Pauli blocking but mean- field contributions U a (r,t) to 
the energy, i.e., <j) -> 0(0 - ±(k +p) 2 - U a - U b , . . .)) 
all A's become nonzero. The space displacements are 
the same as in (iii). The energy and the momentum 
corrections 2Ae,k — —A t -^(U a + Ub), represent the 
energy and the momentum which the effective molecule 
gains during its short life time At. These corrections 
are in fact of three-particle nature, however, only on the 
mean-field level, thus formally within the binary pro- 
cess. The three scalars, At, d and a, are still sufficient 
to parameterize all A's. We note that in the energy- 
conserving 8 function all A's and mean-fields compensate 

5(El+£2 _ £3 _ £4 _ 2A£) = j(| + £_M!_M!' 



This limit describes the dilute quantum gases and (0) re- 
duces to the kinetic equation found in . 

(v) With the medium effects on collisions, all A's become 
independent. A generally nonparabolic momentum- 
dependency of the quasiparticle energy does not allow to 
separate the center-of-mass motion, and an anisotropic, 
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inhomogeneous and time-dependent distribution in the 
Pauli blocking ruins all symmetries of the collision. The 
energy conservation does not reduce to the simple con- 
servation of the kinetic energy. The uncompensated 
residuum of the energy gain contributes to a conversion 
between the kinetic and the configuration energies of the 
system. 

The kinetic equation (^) is numerically tractable by re- 
cent Monte Carlo codes. Kortemayer, Daffin and Bauer 
pj have already studied the kinetic equation with an in- 
tuitive extension by Enskog-type displacement, claiming 
only a little increase in numerical demands compared to 
the BE. The larger set of A's in (Q) does not require 
principal changes of the numerical method used in 
The functions e, z, |T S ^| 2 , and A's should form a consis- 
tent set, i.e., it is preferable to obtain them from Green's 
function studies, e.g., like Eventual fitting param- 
eters should enter directly the effective nucleon-nucleon 
interaction. 



IV. OBSERVABLES 

Let us presume that a reasonable set of functions e, 
z, |T S ^-| 2 , and A's is known. One can then proceed to 
evaluate observables. Here we present the density n a of 
particles a, the density of energy £, and the stress tensor 

Jij- 

The observables in question are directly obtained from 
balance equations which also establish conservation laws. 
Integrating the kinetic equation (j^) over momentum k 
with factors e±,k, 1, one finds that each observable has 
the standard quasiparticle part following from the drift 
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and the A-contribution following from the nonlocality of 
the scattering integral 
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where P = |T«| 2 2^( £l +£ 2 -e 3 -e4)/i/2(l-/3-/4). The 
arguments denoted by numerical subscripts are identical 
to those used in ([j]), for all A's are explicit. The T-matrix 
T cx used in the quasiparticle part of energy is the real 
part of the antisymmetrized Bethe-Goldstone T-matrix, 
= (1 - 6 ab )T£ + 6 ab V2T£. The actual observables 
are sum of the quasiparticle part and the A-correction. 

In the low density limit, the A-contributions (||) be- 
come proportional to the square of density. Therefore 
they turn into the second order virial corrections. In the 
degenerated system, the density dependence cannot be 
expressed in the power-law expansion since the T-matrix, 
and consequently all A's, depend on the density. Never- 
theless, we find it instructive to call the A-contribution 
the virial corrections because of their similar structure. 

The total energy £ — £ qp + AS conserves within ki- 
netic equation (||) . This energy conservation law general- 
izes the result of Bornath, Kremp, Kraeft and Schlanges 
fl9f restricted to non-degenerated systems. At degener- 
ated systems, a new mechanism of energy conversion ap- 
pears due to the medium effect on binary collisions. This 
mechanism can be seen writing the energy conservation 
" £ 



W £qP + W A£ = in a form 



d£ qp 
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The first term on the right hand side is the drift contri- 
bution to the energy balance. It is the only mechanism 
which appears in the absence of the virial corrections. 
The second term is the mean energy gain. It provides 
the conversion of the interaction energy controlled by 
two-particle correlations into energies of single-particle 
excitations. By this mechanism, the latent heat hidden in 
the interaction energy is converted into "thermal" excita- 
tions. The non-zero energy conversion, Ae ^ 0, results 
from the time-dependency of the scattering phase shift 
on the quasiparticle distribution via the Pauli blocking 
of internal states. Accordingly, the energy conversion is 
a consequence of the in-medium effects. 

The energy gain has its space counterpart in the mo- 
mentum gain to the stress forces -SpJi^. The energy 
density contributing to the stress tensor (Q) has the gra- 
dient 

dkdpdq 



d£ qp _\r^ f dk dfi 

a K ' a,b 
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In the absence of the virial corrections (B), the first term 
of ([n]) combines together with the derivative of the first 
term in (|^) into the standard quasiparticle contribution. 
In the presence of the virial corrections, the energy gain is 
necessary to obtain the correct momentum conservation 
law. 

The density of energy £ = £ qp + A£ given by (||) and 
(^) alternatively results from Kadanoff and Baym for- 
mula, 
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2^2 



*1 
2m 
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with G < in the extended quasiparticle approximation (0). 
The particle density n a — n^ p + An a obtained from (5) 
via the definition, n a = J |^ j0js G < i a l so confirms (£) 
and (||). The equivalence of these two alternative ap- 
proaches confirms that the extended quasiparticle ap- 
proximation is thermodynamically consistent with the 
nonlocal corrections to the scattering integral. 

For equilibrium distributions, formulas (ph and (0) pro- 
vide equations of state. Two known cases are worthy of 
comparison. First, the particle density n a = + An a 
is identical to the quantum Beth-Uhlenbeck equation of 
state [p^pil , where n^ p is called the free density and 
An a the correlated density. Second, the virial correction 
to the stress tensor has a form of the collision flux con- 
tribution known in the theory of moderately dense gases 



V. SUMMARY 

In this Letter we have derived the kinetic equation (7) 
which is consistent with thermodynamic observables (£) 
and up to the second order virial coefficient. The 
presented theory extends the theory of quantum gases 
p3| , pT[ and non-ideal plasma Q to degenerated system. 
The most important new mechanism is the energy con- 
version which follows from the medium effect on binary 
collisions. 

The proposed corrections can be evaluated from known 
in-medium T-matrices and incorporated into existing 
Monte Carlo simulation codes, e.g., with the routine used 
in ||. The nonlocal corrections to the scattering inte- 
gral and corresponding second order virial corrections to 
thermodynamic observables enlighten the link between 
the kinetic equation approach and the hydrodynamical 
modeling. With expected progress in fits of the single- 
particle energy and other ingredients of the kinetic equa- 
tion, the virial corrections can improve our ability to infer 
the equation of state from the heavy ion collision data. 
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